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Abstract 

In general relativity, the generators of symmetries are called Killing vector fields. The spacetime 
geometry of a system of two point masses moving on a circular orbit has a helical Killing vector 
(HKV). We show how Kepler's third law for circular orbits, and its generalization in post-Newtonian 
theory, can be recovered from a simple, covariant condition on the norm of that HKV. This unusual 
derivation can be used to illustrate some concepts of prime importance in a general relativity course, 
including those of Killing vector, covariance, coordinate dependence, and gravitational redshift. 
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I. INTRODUCTION AND SUMMARY 



One could hardly overstate the role played by symmetries in physics. Symmetry consid- 
erations can very often drastically simplify the process of solving a given physics problem. 
At a more fundamental level, symmetries are deeply connected to the existence of conserved 
quantities {via Noether's theorem), and they considerably restrain the span of admissible 
field theories in modern theoretical physics. Within Einstein's theory of general relativity 
(GR), which describes gravity as a manifestation of the curvature of spacetime, the gener- 
ators of symmetries are called Killing vectors. They are widely used in current research in 
gravitation theory, and are an essential part of one's education in GR.^"^ 

Gravitational radiation is one of the most important predictions of GR. Observing gravi- 
tational waves would have a tremendous impact on physics, astrophysics and cosmology.^ A 
worldwide effort is currently underway to achieve the first direct detection, using kilometer- 
scale ground-based interferometers such as LIGO,^ Virgo, ^ GEO600,^° and LCGT,^^ as well 
as future space-based antennas, such as the planned ESA-led mission eLISA/NGO.^^ Binary 
systems composed of compact objects (neutron stars or black holes) are among the most 
promising sources of gravitational radiation. However, the detection and analysis of these 
exceedingly weak signals require very accurate theoretical predictions, for use as template 
waveforms to be compared to the output of the detectors. Hence, modeling the orbital 
dynamics and gravitational-wave emission of compact binary systems is a timely problem 
in relativistic astrophysics. 

Except for the occurrence of a gradual inspiral driven by gravitational radiation-reaction, 
the orbits of stellar-mass compact binaries can be considered to be circular, to a very high 
degree of approximation.^^ Mathematically, the approximation of an exactly closed circu- 
lar orbit translates into the existence of a helical Killing vector K'^, along the orbits of 
which the spacetime geometry is invariant. Given the high astrophysical relevance of this 
approximation, helically symmetric spacetimes have been studied extensively. ^^"^^ 

In this paper, we shall consider two non-spinning compact objects moving on exactly 
circular orbits. These will be modeled as point masses rriA (with A = 1,2), a prescription 
commonly adopted in the field of gravitational- wave source modeling. '^^''^^ We will prove 
that the gradient of the norm of the helical Killing vector K"' must vanish along the 
worldlines of the particles: 



We will then show how this simple, geometric result can be used to derive the main relation 
encoding the orbital dynamics of the binary, namely Kepler's third law (for circular orbits), 
and its generalization in post-Newtonian theory. This unusual derivation can be used to 
illustrate numerous concepts of prime importance in a GR course, including those of Killing 
vector, covariance, coordinate dependence, and gravitational redshift. 

This paper is organized as follows: In Sec. II we summarize some well-known, yet useful 
properties of Killing vectors. We then discuss, in Sec. Ill, the physics of binary systems 
of point masses moving along circular orbits, introducing the notion of redshift observable 
in Sec. Ill A, and proving the relation (1.1) in Sec. IIIB. The derivation of Kepler's third 
law, and its generalization in post-Newtonian theory, are discussed in Sees. IV A and IV B, 
respectively. Finally, Sec. V is devoted to some further comments of physical relevance on 
the norm of the HKV and its link to the redshift observable. 
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II. SOME PROPERTIES OF KILLING VECTORS 



For the convenience of the reader, we start by summarizing a few well-known, elementary 
properties of Killing vectors, which will be used extensively throughout this paper. The 
main property of a Killing vector field fc" is that it satisfies Killing's equation 

£kgal3 = '^aki3 + Vf3K = 0, (2.1) 

where Vq is the covariant derivative compatible with the spacetime metric ga/3, and £k is the 
Lie derivative along fc". Equation (2.1) expresses the invariance of the spacetime geometry 
along the integral curves of the Killing vector. 

In addition. Killing vectors straightforwardly provide well-defined conserved quantities: 
Let be the four-velocity of a test particle, tangent to its worldline, and normalized such 
that Qapu^u" = —1. Then the scalar product s = k" constant of the motion along 

the timelike geodesic followed by the test mass: 

s = u^Vp {u"k^) = ii'^k^ + wVV/jfc^ = . (2.2) 

This immediately follows from the geodesic equations of motion, = 0, and the antisym- 
metry property of the tensor Vq/c/j. Two familiar examples of such conserved quantities in 
curved spacetime are the energy per unit mass e = — /c^-jMq, and angular momentum per unit 
mass j = k'^^^Ua of a test particle in orbit around a rotating black hole, where fc^^ and /c^^ 
are Killing vectors associated with the stationarity and axisymmetry of the Kerr metric. 



III. HELICALLY SYMMETRIC POINT-PARTICLE SPACETIMES 

We now consider a binary system of non-spinning compact objects moving on a circular 
orbit. The neutron stars or black holes will be modeled as point particles with constant 
masses (with A = 1, 2), and four- velocities u'\ normalized to = — 1- Note that 

the point masses tua are not test particles; their stress-energy tensor curves the geometry 
via Einstein's field equations. Their motion obeys the standard geodesic equations, u\ = 0, 
albeit in a regularized metric g^^ such that the divergent self-fields of the point particles 
have been carefully subtracted. ^^"^^ 

The spacetime geometry of that binary system is neither stationnary, nor axisymmetric; 
however it is invariant along the orbits of a helical Killing vector (HKV) i^T". Far away from 
the binary system, this field has the asymptotic behavior 

K--^{d,r + Q{d^f, (3.1) 

where the four-vectors (9j)" and (<9^)" are part of the coordinate basis of an inertial frame of 
reference. The constant Q is interpreted as the circular-orbit frequency of the binary system. 
Heuristically, K°' can be seen as the generator of time translations in a co-rotating frame. 
In particular, if you imagine yourself "sitting" on one of the particles, orbiting around the 
companion star, then you would observe no change in the local geometry. In other words, the 
metric is invariant along the worldlines of the particles: £uA9ap = 0. This implies that the 
four-velocities of the particles must be aligned with the HKV evaluated at their respective 
coordinate locations: 

u'X = ulK'X ^ K^ = ZAul, (3.2) 
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FIG. 1. Spacctime diagram picturing a binary system of point masses ruA — 1,2) on a circular orbit 
with constant azimuthal frequency fl. The hehcal KiUing vector K" — (9t)" + {d^)'^ is ahgned with the 
four-velocities tangent to the worldlines Ta of the particles. The redshift observables are given by the 
scalar products za — —K'^u^. 

where we pose za = See Fig. 1 for an illustration. As a coefficient of proportionality 

between two four-vectors, must be a scalar. It can be assigned several physical inter- 
pretations, which, in a GR course, would provide the opportunity to discuss the important 
notion of coordinate invariance. 



A. Redshift observable 

First, contracting the relation (3.2) with the four-velocity u^, and remembering Eq. (2.2), 
we notice that za is a constant of the motion associated with the helical symmetry: 

Za = —KaW^ = const. (3.3) 

Furthermore, it can easily be established that za measures the redshift of lights rays emitted 
from particle A, and received far away from the binary, along the helical symmetry axis 
perpendicular to the orbital plane:^^ Let be a four-vector tangent to the worldline of 
such a light ray, e.g. the four-momentum of the associated "photon". Then the ratio of the 
photon energy at reception and emission is given by 

Srec ^ (W>«)rec _ (^>a)rec _ ^ 

We made use of the equality w^g^ = K^^^ = {dtY between the four-velocity of the observer 
and the HKV, of the relationship (3.2) at the location of the emitter, and of the conservation 
of K°'Pa along the null geodesic of the photon. (The proof of that last point is identical to 
that given in Eq. (2.2), with the substitutions — )■ and fc" — )■ K'^.) See Fig. 2 for an 
illustration of this Gedankenexperiment. Following Detweiler, we shall thus refer to za as 
the "redshift observable". 
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FIG. 2. A photon 7 with four- momentum is emitted from the particle mi with four- velocity u"^ = uf, 
and received far away from the binary system, by a distant inertial observer with four- velocity m"(,c = i^t)"', 
along the helical symmetry axis z perpendicular to the orbital plane. Detweiler's observable zi measures 
the redshift £rcc/^^cm of the photon, or equivalently the ratio dri/dt. 



Finally, in a cylindrical coordinate system {ct, p, z} adapted to the helical symmetry, 
i.e. such that the expression K"" = (9^)" -|- Q (9<^)" holds everywhere, and not merely far 
away from the binary, Eq. (3.2) implies 

z^=K)-^ = ^. (3.5) 



Hence the redshift observable coincides with the inverse time component of the four-velocity 
of the particle, or equivalently with the ratio of the proper times elapsed along the worldlines 
of the particle and of the distant inertial observer (c/. Fig. 2). This last interpretation is in 
agreement with the usual notion of redshift (gravitational redshift and/or Doppler effect). 



B. Geometric characterization of the binary dynamics 



We now have at our disposal all of the concepts and results necessary to prove Eq. (1.1), 
namely that the spacetime gradient of the norm = g^pK^^K^ of the HKV must vanish 
at the location of each particle. The derivation goes as foUows:^^ 

^ (V„i^2)^ = Rf'^iV^Kp)^ = -ZAU^A^pK^)^ = -ZA {zAut + ZAiit) = . (3.6) 

We successively made use of the relationship (3.2) between K'^ and u% of Killing's equation, 
of the geodesic equations of motion, and of the fact that za is a constant of the motion. 
Although pretty straightforward, that proof would provide a good exercice for students. 

The formula (1.1) is quite simple and elegant; it is covariant, and only makes reference to 
well-defined geometrical concepts in GR. It implies, in particular, that the function -ft'^(x) 
has extrema at the coordinate locations of the particles (see Fig. 3). Furthermore, we are 
going to show that (1.1) encodes a well-known result of classical mechanics, namely Kepler's 
third law (for circular orbits), and its generalization in post- Newtonian theory. 
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IV. KEPLER'S THIRD LAW AND ITS RELATIVISTIC GENERALIZATION 



A. Newtonian gravity 

For simplicity, we start by restricting ourselves to the Newtonian approximation of the full 
theory of GR. This corresponds to the leading-order results in the formal limit c"^ 0. In 
that weak-field, non-relativistic approximation, the spacetime metric expressed in cylindrical 
coordinates {ct, p, ip, z} takes the form 



ds^ = ( -1 + ^ ) cMt' + dp' + p'dip' + dz' , (4.1) 



where the Newtonian gravitational potential reads t/(t,x) = ^^Gm^/jx — y^it)], with 
yA{t) the coordinate trajectory of the mass rriA- In our convenient coordinate system such 
that = {dtf + n {d^f, the norm of the HKV reads K' = goo + 2fi^o^/c + n'g^^/c\ 
With the explicit expression (4.1) of the metric, this gives 

K' = -l + ^ + ^ . 4.2 

In order to apply the result (1.1), we need first to compute the partial derivatives of the 
scalar K'. Using the Cartesian coordinates {x*} associated with the cylindrical coordinates 
{p,(p, z} in the usual way, namely = pcosyj, = psimp, and = z, we have 

dtK' = ^dtU and d,K' = ^{d,U + pQ'n'), (4.3) 

where n = (cos<^, sin<^, 0) is the unit vector in the orbital plane z = 0. Focusing first on the 
spatial components, the relationship (1.1) implies 

-p^i]v;, = (9,t/)^. (4.4) 

We recognize Newton's third law expressing the equality of the centripetal acceleration of 
body A and of the Newtonian gravitational force exerted on the body A by the body B ^ A. 
Computing the forces explicitely, we find^® 



Pin' = ^ and p^n' = —^, (4.5) 



where r = |yi — y2| = pi + p2 is the coordinate separation between the two point masses. 
We may then add up Eqs. (4.5), or remember that in the center-of-mass frame pi = rm2/m 
and p2 = rrrii/m, with m = mi + m2 the total mass of the binary. Both solutions yield 
Kepler's third law 

tf-^. («) 

which is recovered here in the particular case of circular motion. (Our derivation cannot be 
extended to generic eccentric orbits, for which the helical symmetry is lost.) On the other 
hand, the time component of Eq. (1.1) is identically satisfied, as {dtU)j^ oc ■ ni2 vanishes 
for circular orbits, = dys/dt being the coordinate velocity of the particle B ^ A, and 
ni2 the unit vector pointing from 1712 to mi. 
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B. Post-Newtonian gravity 



Equation (1.1) being valid beyond the Newtonian limit, the result (4.6) can be extended 
to include corrections coming from the full theory, by keeping known post-Newtonian terms 
in the metric Qafs- For example, at the first post-Newtonian (IPN) order, i.e. including the 
general relativistic corrections 0{c~'^) to the Newtonian expression (4.1), the metric reads 
(in Cartesian-like harmonic coordinates)^^ 



fl'oo 



2Gmi 1 
-1 + ^— + ^ 



-|- G'^mim2 



— (4vi ■ vi - (m ■ w,f) - 
'1 '1 



rir2 



^' 12 



+ 



2r2rf2 2r2ri2/_ 



+ (1^2) + 0(c-6) 



9i3 



AGmi 



iGm^ 



6'^ 1 



2Gmi ^ 2Gm2 



+ 0{c 



(4.7a) 
(4.7b) 
(4.7c) 



where 6^^ is the usual Kronecker symbol, ri2 = |yi — y2| = t is the coordinate separation, 
ta = |x — y^l, and = (x — yA)/fA- Transforming that metric to cylindrical coordinates, 
and repeating the calculation detailed in Sec. IV A, we recover from {diK'^)y^ = the known 
generalization of Kepler's third law at IPN order (in harmonic coordinates), namely^^ 



Gm 



1 + (-3 + z.) ^ + Oic-'] 



(4. 



At that order of approximation, the calculation involves a crucial contribution coming from 
the general relativistic frame-dragging effect, through the metric component go^p- The IPN 
coefficient (—3 + u) in Eq. (4.8) depends on the symmetric mass ratio u = mim2/m? , such 
that 1/ = 1/4 for an equal-mass binary, and z/ — )• in the extreme mass ratio limit. (The 
relation {dtK'^)j^ = is still found to be satisfied identically.) 

In a GR course, this derivation would provide the occasion to discuss the key notions of 
covariance and coordinate dependence using a concrete example: Being covariant, Eq. (1.1) 
conveys a physically meaningful result, independent of a particular choice of coordinates. By 
contrast, the generalized version (4.8) of Kepler's third law is coordinate dependent; the IPN 
coefficient could be different from (—3-1- 1^) if the relationship between the invariant frequency 
VL and the coordinate dependent separation r was expressed in another coordinate system. 
But the precise way in which the function VL{r) changes, depending on the coordinate system 
used to write the post-Newtonian metric, is precisely encoded in the covariance of Eq. (1.1). 
Finally, we note that while the relationship between VL and r is coordinate dependent, the 
functions za{^) are coordinate invariant; they have recently been computed up to very high 
orders in the post-Newtonian approximation. ^^"^^ 



V. HELICAL KILLING VECTOR AND REDSHIFT 

From the Newtonian result (4.2), we have the asymptotic behavior ~ {pVL/cY > 
in the limit p — )■ -|-oo, which indicates that the HKV is spacelike far away from the helical 
symmetry axis. Close to the binary system however, and along the worldlines of the particles 
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FIG. 3. The norm of the hehcal Kilhng vector as a function of the coordinate w along the direction 
joining the two particles, within the orbital plane. 



in particular, the HKV is timelike {K'^ < 0) [remember Eq. (3.2)]. In between, there must 
exist a worldtube over which the HKV is null {K"^ = 0). This hypersurface is usually referred 
to as the "light cylinder". In the flat spacetime limit rriA — )■ 0, its radius is simply p = c/fl. 
(In the langage of classical mechanics, this is the distance from the axis of rotation for which 
the velocity of an observer rotating at the constant angular rate fl reaches the vaccum speed 
of light c.) 

Furthermore, when evaluated at the locations of the particles themselves, the norm of the 
HKV is directly related to the redshift observables. Indeed, from Eq. (3.2) we immediately 
get z\ = —K\. In the Newtonian limit, we thus have 

(5.1) 

This result is consistent with the interpretation of za as a measure of the redshift of light rays, 
as discussed in Sec. HI A. The observable za has two contributions: (i) a term proportional 
to the Newtonian potential Ua evaluated at the coordinate location of particle A, which gives 
the gravitational redshift, or Einstein effect, and (ii) a term involving the relative velocity 
va = Pa^ with respect to the distant observer, yielding a transverse Doppler effect. (In the 
flat spacetime limit niA — 0, we recover the special relativistic result za = a/1 — v\/c'^.) 

Based on the previous quahtative discussion, the function K'^{w) is depicted schematically 
in Fig. 3, with w = pcos {ip — pi) the coordinate along the direction joining the two particles, 
within the orbital plane. (The divergent Newtonian self-fields of the point masses are shown 
in dashed lines; these are well-known artifacts of the use of point particles to model the 
actual physical compact stars, which are extended objects.) A more quantitative analysis 
of the function K'^i'x) in Newtonian (or post-Newtonian) gravity, in and out of the orbital 
plane, could be a useful exercise for students. In particular, they could be asked to plot that 
relation for different sets of values for {fl, mi, 1712}- 
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